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Abstract. We present a topological proof of the existence of a normally hyperbolic 
invariant manifold for maps. In our approach we do not require that the map is a 
perturbation of some other map for which we already have an invariant manifold. 
But a non-rigorous, good enough, guess is necessary. The required assumptions are 
formulated in a way which allows for rigorous computer assisted verification. We apply 
our method for a driven logistic map, for which non-rigorous numerical simulation in 
plain double precision suggests the existence of a chaotic attractor. We prove that this 
numerical evidence is false and that the attractor is a normally hyperbolic invariant 
curve. 
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1. Introduction 

In this paper we give a proof of existence of normally hyperbolic invariant manifolds for 
maps. The construction is performed in the state space of the map. Assumptions needed 
for the proof are of twofold nature. First we require topological conditions which follow 
from suitable alignment of the coordinates (these are the so called covering relations). 
Next we require that our map satisfies cone conditions. The aim of the paper though is 
not to produce yet another proof of the normally hyperbolic invariant manifold theorem. 
Our aim is to produce a tool that can be applied in rigorous-computer-assisted proofs. 
To show the strength of our approach we apply our theorem to a driven logistic map 
introduced in [2] . The considered map is such that standard numerical simulation gives 
evidence of a chaotic attractor. The example is a demonstration of the fact that one 
has to be careful with the arithmetics in simulations, since the numerical evidence of 
an attractor is false. The map in fact possesses a normally hyperbolic invariant curve. 
This is apparent when simulations are performed using multiple precision computations. 
The strength of our method lies in the fact that even for such an example, which defeats 
standard numerical simulations, we are able to produce a rigorous proof of existence of 
a normally hyperbolic invariant curve. 

The approach to normally hyperbolic manifolds presented here is in the spirit of [3] 
and [6]. In ^ a topological proof of existence of invariant sets with normally hyperbolic 
type properties is given. In [6] the result is extended to prove normally hyperbolic 
invariant manifolds. In both cases the proofs relied on assumptions that the first iterate 
of the map is well aligned with the stable and unstable manifolds. Similar approach 
was also used in j^J to give a proof of existence of a center manifold. The result in 
[1] is for ODEs and relies also on the fact that hyperbolic dynamics is uniform. The 
main difference between our paper and results mentioned above is that we assume that 
hyperbolic expansion and contraction aligns with the tangent spaces of the invariant 
manifolds after a suitable (possibly large) number of iterates of the map. This setting 
is more general, and also more typical for normal hyperbolicity. 

The paper is organized as follows. Section [2] introduces basic notations used 
throughout the paper and provides a setup and an outline of our problem. Section 
E] contains a geometric construction of a normally hyperbolic manifold. We first give 
a construction of a "center-stable" manifold (the term "center-stable" refers to the 
normally hyperbolic invariant manifold union its associated stable manifold; analogous 
terminology is used by us for the "center-unstable" manifold). A center- unstable 
manifold is obtained using a mirror construction to the center-unstable manifold, by 
considering the inverse map. The intersection of the center-stable and center-unstable 
manifolds gives us the normally hyperbolic invariant manifold. In Section H] we show 
how to verify assumptions of our theorems using local bounds on derivatives of the map. 
In Section [5] we present our example of the driven logistic map and apply our method 
to it. 
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2. Setup 

We start by writing out some basic notations which we shall use throughout the paper. 
A notation Bi{q,r) will stand for a ball of radius r centered at q in W. We will also 
use a notation Bi — 5^(0, 1). For a set A we will denote by A its closure, by int^l its 
interior and by dA its boundary. For a function / we will use a notation dom(/) to 
denote its domain. For points p — {x,y) we shall use notation Trx{p), T^y{p) to denote 
the projection onto the x and y coordinates respectively. 

We now introduce the setup of our problem. Let D and U be open subsets in 
such that D cU. Let 

f:U^U, 

be a diffeomorphism. Let m, s, c e N be such that u + s + c = n. We assume that there 
exist a diffeomorphism 

: W ^ C M" X X A 

such that (f){c\D) = := Bu x Bs x A, and A is a compact c dimensional manifold 
without boundary. We define f^:D^^WxWxAas 

f<t>^(t>° f ° 

We assume that there exists a finite covering {Ui}i^i of A and an atlas 

r]i:Ui^ Be- 
Throughout the work wc will use a notation 

B^B^xB.x B,. 

For we consider local maps fji : B D dom(/y) — >■ x x defined as 

fij — Vj °f4>° 

fji := (id, id, r)i) for i e I. 

Note that the domain of fij can be empty, and will usually be smaller than B. The 
following graph depicts the above defined functions and their mutual relations. 

D 

Bu X Bg X A 
B 

Our task in this paper will be to find a normally hyperbolic invariant manifold, 
together with its stable and unstable manifolds within the set D. 

We will use the following notations for our coordinates: x e R", y G W, 9 G Be, 
A G A. The coordinate x will play the role of a globally unstable direction, and the 
coordinate y will play the role of a stable direction for the map (hence the superscripts 



fji 



u 

X R" X A 

r X R" X 'Br 
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u and s, which stand for "unstable" and "stable" respectively). The coordinate A will 
play the role of the central direction, in which the global dynamics is weaker than in 
the stable and unstable coordinates. The notation 6 will also be used for the central 
direction, but it will be reserved to denote the central coordinate in the local coordinates; 
i.e. 9 = rii{X) for some A G A and i E I. 

3. Geometric approach to invariant manifolds 

In this section we give the construction of a normally hyperbolic invariant manifold. The 
construction is performed in the state space of our map. It is based on the assumptions 
of covering relations and cone conditions. We first give an introduction to these tools in 
Section [STTl In Section [3^ we formulate our assumptions on the map in terms of covering 
relations and cone conditions, which will imply the existence of a normally hyperbolic 
manifold. In Section 13.31 we show how to construct a center-stable manifold of our 
map. The construction of a center-unstable manifold follows from a mirror argument. 
The intersection of center-stable and center-unstable manifolds gives us a C° normally 
hyperbolic invariant manifold. Let us write explicitly that for a normally hyperbolic 
manifold which does not have an associated stable manifold, the center-stable manifold 
will be the the normally hyperbolic manifold itself. Analogous statement holds also for 
center-unstable manifolds. 

3.1. Covering relations and cones 

Covering relations are topological tools used for proofs of nontrivial symbolic dynamics 
of dynamical systems. The method is based on the Brouwer fixed point index, and 
the setting is such that it allows for rigorous numerical verification. The method has 
been applied in computer assisted proofs for the Henon map, Rossler equations [19], 
[6], Lorenz equations [8], Chua circuit [7] or Kuramoto-Shivashinsky ODE [18], amongst 
others. The method is based on singling out a number of regions, called h-sets, which 
have hyperbolic type properties. Using these properties one can find orbits of the system, 
which shadow the h-sets along their trajectories. The method of covering relations 
relies on the system having expanding and contracting coordinates. In this section we 
generalize covering relations to include also a central direction. The setup is similar to 
that of [3], 0, but has been simplified. Our proofs are now simpler and based only on 
continuity arguments. They no longer require the use of degree theory, with little loss 
of generality. 

For any p = (x, 6*) G B and r„, r^, > we introduce a notation 
N{p,ru,rs,rc) := B^.(x,r„) x Bs{y,rs) x Bc{9,rc). 

We define 

N- = N~{p,ru,rs,rc) := dBu{x,r^) x Bs{y,rs) x Bc{e,r^) 
= N'^{p,ru,rs,rc) 

:= 5„(x,r„) X ((R^ x \ {B,{y,r,) x 5,(0, r,))) . 
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We assume that all boxes N which we are going to consider here are contained in B. 
We will refer to a box as a ch-set (center-hyperbolic set) centered at p. 

In following arguments we shall often consider different ch-sets. To keep better 
track of our notations and to make our arguments more transparent we shall stick to 
a convention that for two ch-sets Ni,N2 centered respectively at pi = {xi,yi,9i) and 
P2 = {x2,y2,02) we shall write 

N, = N,ip,,rlrlri) := x BUy^rl^ x XiOi,ri) for z = 1, 2. 

Definition 1. Let (7 : B — t- R" x x 5^ be a continuous function. Let pi = (xj, yi, 9i) 
fori = 1,2 and let Ni, N2 be two ch-sets in B centered at pi and p2 respectively. We 
say that Ni (7-covers N2 if 

g{p,) e int(iV2), (1) 
nMNn)^Blix2,rl) = (D, (2) 
g{N,)nN+ = iD. (3) 

In such case we shall write Ni =^ N2. 

Remark 2. Definition [1] is a simplified definition of a covering relation. More general 
versions can be found in [8], [9], [19] in the setting of hyperbolicity, or in [3], [5] in a 
setting when additionally a central direction is included. 

For 7 = (a, 6, c) e R^, and q = (x, y, 6*) G R" x R"* x R^ we define 
: R" X R" X R'^ ^ R 

Q^{q) ■.= a\\xf + b\\yf + c\\9f . (4) 

If a > 6, c < 0, then for p e R" x R"* x R'^ we will refer to 

C{p,l):={q:Q,{p-q)>0} 
ClS db horizontal cone centered at p (see Figure [2]). 

Definition 3. Let N be a ch-set and 7 = (a, b, c) be such that a > 0, b,c < 0. We will 
refer to a pair (A^, 7) as a ch-set with cones. 



Computer Assisted Proof for Normally Hyperbolic Invariant Manifolds 6 

Definition 4. Let (A^, 7) = {N{{x,y,6),ru,rs,rc),'~f) be a ch-set with cones. A 
continuous function h : Bu{x,ru) N is called a horizontal disc in (A^, 7), iff 
7r^h(x) = X and for any x*,x** G Bu{x, r„), 

Q^{h{x*)-h{x**))>0, (5) 

Lemma 5. Let Ni = Ni{{xi,yi,9i),rl,rl,rl) fori = 1,2 and let (A^i,7i), (A^2,72) be 
two ch-sets with cones. Assume that 

N^ ^ N2 (6) 

and that for any q*, q** G A^i such that q* ^ q** and Q^^{q* - q**) > we have 

Q,M<il-9{q*l)>^- (7) 

//hi is a horizontal disc in (A''i,7i) then there exists a horizontal disc h.2 in {N2,'~f2) 
such that g{h,(Blixi, rl))) n N2 = h2(5'(x2, r^)). 

Proof. Without loss of generality we assume that Pi = P2 = and that = 1 for 
i = 1,2 and n G {u, s, c}. In other words we assume that for i = 1,2 

AT, = si X ^1 X si = S„(0, 1) X Bs{0, 1) X 5,(0, 1). 

Let 7j = (oj, hi, Cj) for i = 1,2 and let h be any horizontal disc in Ni. Then by dl]), 
(ED and dZD for 2;*,a;** G bI, x* 7^ x** 

a2 WTr^Hx*)) - 7r.(?(h(x**))f > g,,(^7(h(x*)) - giHx**))) > 0, (8) 

which means that tt^ o g o h is a. monomorphism. 

Using a notation hi(x) = (x, hi{x)) G 5^ x (5* x 5^), for a G [0, 1], we define a 
family of horizontal discs hQ,(x) = {x, ahi{x)). Let Fq, : B^ — )■ be a continuous family 
of functions defined as 

F„(x) ■.= n^ogo ha{x). 

We shall show that B^ C Fi(i?^). Functions Fq, are monomorphisms, hence sets 
Ac := Fa{BD are homeomorphic to balls in M"; moreover c^Aq, = Fa{dBD. By 
Definition m of a horizontal disc, h^idBl) C A^^f. From assumption by conditions 

©, (ED 

vrx.^7(0) G S^, (9) 

n 5' c F,(A^r) n fi' = 0. (10) 

From the fact that G i?,', 

Fo(0) G Fo(S^) = Ao. (11) 
Since ho(0) = 0, by (P 

Fo(0) = n,ogo ho(0) = 7r,(7(0) G bI, (12) 

From (dlD, ([nD follows that AqHbI^ 0. This by (HDD implies that 5^ C Aq. By 

2 

continuity of with respect to a this means that B,^ C A„ for all a G [0, 1]. In 
particular B^ C Ai = Fi{Bl). 
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Figure 2. Covering relations for two iterates of a map /. For the second iterate of 
the map the coordinate x is expanding and y is contracting (for the first iterate of / 
they are not). The fact that expansion in x is stronger than expansion in is visible 
from the fact that the cones C{f'^{p),^3) are "tighter" than cones C{p,ji). 



2 2 

Since Fi is a monomorphism and C Fi{BD, for any v E there exists a unique 
x{v) G Bl such that Fi{x) = v. We define h.2{v) = {v, h2{y)) := (f , 'Kyfiogoh.i{x{v))). 



For any v* ^ v**, v*,v** G B^, by ([5]) and ([7]) we have 

Q,, (h^iv*) - h2{v*)) = Q,,{g o h,{x{v*)) - 9 ° M^iv**))) 

>g^,(hi(x(t;*))-hi(x(t;**))) 
> 0. 

Since Q^, (h2(t;*) - h2(t;**)) > 

02 ||f * — V** II 

> -&2 Ik, (h^iv*) - h^iv**))^ - C2 llvr, {h2{v*) - h2{v**))f 

> min(-62, -C2) \\h2{v*)-h2{v**)\\\ 
and therefore h2 is continuous. 



□ 



Remark 6. Let us note that since we have freedom of choice of the radii r„, and Tc 
it is not necessary for x to be expanding, y to be contracting and 6 to have weaker 
dynamics for each single iterate of the map. In Figure |2] we have a sketch of a situation 
in which x becomes expanding and y contracting after a second iterate. In Figure |2] the 
coordinate 6 is expanding. It will turn out that such a scenario is acceptable for us and 
can be dealt with by increasing Tc for successive iterates. 



3.2. Covering relations and cone conditions for normal hyperbolicity 



In this section we formulate our assumptions which will imply the existence of a 
normally hyperbolic manifold. The assumptions are in terms of covering relations 
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and cones and are in the spirit of [5j. There are two major differences though. The 
first is that assumptions used in [5] required the system to have uniform expansion 
and uniform contraction for the first iterate of the map. Here we set up our 
coordinates in the directions of global contraction and global expansion. In the setting 
of normal hyperbohcity the coordinates of global contraction and expansion need not be 
contracting and expanding for the first iterates of the map. What is important is that 
they dominate after a sufficiently large numbers of iterates, in other words, that the 
Lyapunov exponents are negative or positive, respectively. We set up our assumptions 
so that they allow for such setting. The second difference is that our setup has been 
significantly simplified with comparison to [5] . This resulted in a slight loss of generality 
(we do not formulate our assumptions in terms of vector bundles as in [5]) but we need 
to consider fewer assumptions. 

Let 1 > -R > p, r > 0. Assume that there exists a finite sequence of points 
Afc G A, G N such that for any k the set I{k) = {i : Bc{r]i{Xk), p) C Bc{0,R)} is 
not empty. What is more, assume that there exists a set J C {(i, k)\i G I{k)} such that 
i,k)ejVi {Bc{Vi{^k),p))- For points {i,k) G J we define sets 

We will need to assume that the points are sufficiently close to each other. We 
will also need to assume that R and p are sufficiently large in comparison to r. This is 
summarized in Assumption [71 The idea behind it is demonstrated in Figure [3], which 
might provide some intuition. 

Assumption 7. Let m > 1 and let 7o = (^i-O) bo, Cq) G M"', 7^ = (ai, bi, Ci) G M.^ satisfy 
S-m > 0, hjn, Cm < for m = 1,2. Let us also define a set M C B as 

M:=5„(0,r) X S,(0,r) X 5,. (13) 

We assume that for any horizontal disc h in a ch-set with cones (M, 7^) and for any 
i E I there exists {l, k) E J such that h(i?u(0, r)) C dom(r/tO?)r^). In addition we assume 
that for any q*, q** in dom(f/t o f]^-^) such that Qj-^{q* — q**) > we have 

Qfoiv. ° v;\q*) - V. ° vi'iq**)) > ^Q-rM* - r), (14) 

and 

h' := fj^ o T^r^ o h|B^(o,r) is a horizontal disc in {M^^^,^q). (15) 

Assumption [7| ensures that for h in some local coordinates fji we can change to 
coordinates rj^ so that h' := fj^ofj^^ oh lies close to the middle of the set M. Assumption 
[7] is also discussed in Section 14. 3[ where conditions which imply it are given. 
Remark 8. Above we use bold font for 7^ = (a^, bj, Cj), i = 0, 1 to emphasize that these 
are fixed constants, and to distinguish them from other 7 = (a, b, c) in our proofs. 

Definition 9. // for any {i, k) E J there exists a sequence of ch-sets with cones 
(A^i, 7i), . . . , (A'n, 7„) (n can depend on {i,k)) and a sequence zq = i,ii,...,in E I 
such that 

M,,fc =: iVo ^ A^i ^ ^ . . . M, (16) 
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h{Bu{0,r)) 




B,{0,R) 



Figure 3. The change of coordmates o 77"^, a horizontal disc h, and the cones given 
by 7g and 'Ji in different local coordinates. Here, for simplicity, the stable coordinate 
is neglected s = 0. 



\ 



Ni 
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1 X: 1 
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I U- 



-'7 



Figure 4. (see Example [T0|) For the first iterates of the map the ch-sets and cones 
are contracted in the x direction. After a number of steps the expansion in x starts to 
dominate. Note that the coordinate Q is expanding. Since expansion in x is stronger 
than expansion in though, the cones eventually become more flat and their level sets 
Q^^ = c are pulled away from the origin. 



then we say that f satisfies covering conditions. 
// in addition for any qi,q2 € A^^^^i, qi ^ qi, 

Q^i+iifk+ikiQi) - fk+iiM) > Q^Mi - 92) (17) 

for Z = 0, . . . , n — 1, and for 7„ = (a, 6, c) we have 

bi 6 ci c 

ai > a, — > -, — > -, (18) 
ai a ai a 

then we say that f satisfies cone conditions. 

Example 10. This example stands behind the pictures from Figure |H Consider u = c = 1 
and s = 0. Assume that fi^i^ = (Ajj)i,j=i,2 = diag{\,2) fi^i^ = (A|)ij=i,2 = diag{2,l), 
fi^i^ = {Afj)ij=i^2 = diag{5, 2). Let 7q = (1, —1) and = (i, — |). We take ch-sets with 
cones {Ni{{0, 0), r^), 7?), for / = 0, 1, 2, 3 with 

— ^11 ~ ^) 
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B,, 









J 











h(B„(0,r)) 




Figure 5. The sequence of covering relations from Definition |9l together with the sets 
Mi(, -^ti,Kij which are the first step of the inductive construction from the proof 

of Theorem [ni 



70 = 7o, 71 = (45, -16-'), 72 = -i^"'), 73 = (^^^ -^S-'), with 5 = 1 + e. For 
sufficiently small r and e we will have (|T6l) and ( IT71) . For sufficiently small e we also 
have f|T8l) . Assume now that 77^ o r/^^-'^ = diag{l, 1 + This fj^ o r^j^-*^ is taken just as 
a hypothetical example, in order to show that even when a switch to new coordinates 
involves an expansion in the central coordinate the Assumption [7] can easily be satisfied. 
We have 

5 




which means that f|T^ holds for m < 4. 

We now introduce a notation U C for a set 

f/:=5,(0,r) X 5,(0,r) X A. (19) 

The set U will be the region in which we will construct an invariant manifold of points, 
which stay within the set D^p for forward iterations of the map J^. 

3. 3. Existence of a normally hyperbolic manifold - Main result 

In this section we use the assumptions from Section 13.21 to obtain the existence of a 
normally hyperbolic invariant manifold inside of the set U defined in (fT9|) . We start 
with a construction of the center-stable manifold. This is given in Theorem [TTJ The 
existence of an center-unstable manifold follows from mirror arguments for the inverse 
map. The normally hyperbolic manifold is obtained by intersecting the center-stable 
and center-unstable manifolds. This is done in Theorem [131 

Theorem 11. // / satisfies cone conditions then there exists a continuous 
monomorphism V : Bs{0,r) x A ^ U such that 



(i) nyV{y, A) = y, nxV{y, A) = A, 
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(a) for any {y, A) G -Bs(0, r) x A and any n eN 
f;{V{y,X))eD^. 

(Hi) for any q G U such that f^{q) G for all n G N, there exists a (y, A) G -8^(0, r)x A 
such that q = V{y, A), 

(iv) if X*,X** G r]l^{Bc{r]i{\k), p)) for some {i,k) G J then for any y*,y** G Bs{0,r) 
such that {X*,y*) ^ (A**,y**) 

Q-ro ° V{y\ A*) - r/i o V{y*\ A**)) < 0. (20) 

Proof. We take any ?/o ^ -Bs(0, r), Aq G A and (to, /^o) ^ such that Aq G ?7^^"'^(i?c(AK(,, p)) 
and define a horizontal disc ho in M^^^i^^ as 

ho(x) := (x,?/o,r7,o(Ao)). 

Since / satisfies cone conditions, using assumption ( fT6|) and applying inductively Lemma 
Ogives us the existence of indexes zi, . . . , z„j G / and of a horizontal disc hi in (M, 7„J 
such that 

MBu) = {m^, o o r7-/(ho(a;)) G M : x G 5„(0, r), and 
% ° 4 ° e for / = 1, . . . , m}. 

By (HE]) 

Q,,(hi(x*) -hi(x**)) > g,„^(hi(x*) -hi(a;**)) > 0, 

which means that hi is a horizontal disc in (M, 7i). From (fT4ll and (fTSll we know 
that there exists (ti,/€i) G J such that h' := f],^^ o r/r^ o hi is a horizontal disc in 
(Mtj^Ki5 7o)- This in particular means that for fi := fj^^ o /^^ o fj'^, there exists an 
X G 5„(0,r) for which fi(hi(a;)) n M,,,^, ^ 0. By (HZD and ([T1D, for any x* ^ x** such 
that hi(x*),hi(x**) G dom(fi) 

Q,, (fi(hi(x*)) - fi(hi(x**))) > mg^„ (hi(x*) - hi(x**)) > 0. (21) 

Repeating the above procedure inductively (starting the second step with the 
horizontal disc h' and local coordinates given by r/^J we obtain a sequence of points 
Xs G -B„(0,r) and indexes {is, Hs) for s G N such that for 

we have 

f^(hi(xs)) G M,^^^^ ioT w <s. 

Since Bu{0,r) is compact, there exists an xq = xo{yo, Xq) G -B„(0,r) such that 
VZ^ of^(hi(xo)) G U for all s G N. We define V{yo,Xo) := V'^ixoiyo, Xo),yo, Xo). To 
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see that V is properly defined suppose that we have two points Xq 7^ Xg* such that 
fj,^ o f (hi(a;o)), f],^ o f (hi(xo*)) G U for all s G N. Then by ([21]) we obtain 

Q,, (f.(hi(x*)) - f.(hi«))) > mg,„ (f.-i(hi(x*)) - f._i(hi(x**))) 

> . . . (22) 

> m'Q^, (hi(x*) -hi(x**)) 

> 0. 

Since m > 1, (1221) implies in particular that 

IItt^. (fs(hi(xo)) - fs(hi(xo))) II ^ 00 as s ^ 00. 

This is impossible since fs(hi(xQ)) is in M^^ ks? which is a subset of B, which is a 
bounded. 

We now need to show ({20]). Suppose that V{y* , X*),V{y** , \**) G Mi^k and 
<57o (^i ° V{y*y '^*) — fii° V{y**, A**)) > O. Applying estimates analogous to ( j22l) we 
obtain a contradiction. 

Continuity of ^ will follow from the fact that 

Q70 (^.0 ° A*) - fi.o ° ^(2/**, A**)) < 0. (23) 
Since 7g = (ag, bg, Cq) with ao > and bo, Cq < (!23l) gives 
> (r^,„ o Viy*, A*) - r/,„ o V{y**, A**)) 
= ao Wrr^Viy*, A*) - 7r.V(r , A**)f + bo \\y* - y**f 
+ Coh.o(A*)-r^,„(A**)f , 

and therefore 

^o\\7r.V{y*,X*)-7:..V{y**,Xnf 

< min(-bo,-Co) ||(y*,r/.„(A*)) - (y** , //.^(A**))!^ 

□ 

Now we move to proving the existence of the normally hyperbolic invariant 
manifold. First we need a definition. 

Definition 12. We say that f satisfies backward cone conditions if f~^ satisfies cone 
conditions, with reversed roles of x and y coordinates. 

We assume that for / Assumption [7] holds with 7o = "fo'^^- We assume also that for 
/^^ Assumption [7] holds with 70 = ^Jj^^"* (with reversed roles of the x and y coordinates). 

Theorem 13. (Main Theorem) Assume that f satisfies cone conditions for 70°^^" = 
(ao,bo,Co) and backward cone conditions with j^^'^^ = (ao,bo,Co) . // 

|a[,| > |a^| and \hl\ < \h^\ (24) 

then there exist continuous monomorphisms : -Bs(0, r) x A — t- f/, : -8^(0, r) x A — >■ 
U and X : A — )■ [/, such that 

Tiy,xW'{y, A) = (y, A), 7r,,AW^"(l/, A) = (x, A), 7r,x(A) = A, (25) 
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and := x(A) is an invariant manifold for f^, with stable manifold iy*(i?s(0,r) x A) 
and unstable manifold W^{Bu{0,r) x A). 

Proof. Since / satisfies cone conditions, applying Theorem [TT] we obtain W^ly, A) as V. 
Since / satisfies backward cone conditions, once again from Theorem [TT] for f~^ we also 
obtain W^{x, A) as function V. From point |i]in Theorem [11] it follows that ( 125|) holds 
for W and W. 

We shall show that for any A G A the sets W'{Bs{0,r),X) and iy"(5„(0, r). A) 
intersect. Let us define F : 5^(0, r) x Bs{0,r) 5„(0,r) x 5^(0, r) as 

F{x,y) := {n,W'{y,X),nyW\x,X)) . 

Since F is continuous, from the Brouwer fixed point theorem follows that there exists 
an (xo,?/o) such that F{xo,yo) = (a;o,yo) • By (!25|) this means that 

W%yo, A) = {n^W^yo, X),yo, A) = (xq, 7r,l^"(xo, A), A) = l^"(xo. A). 

Now we shall show that for any given A G A there exists only a single point 
of such intersection. Suppose that for some A G A there exist {x* , y*) , {x** , y**) G 
S„(0,r) X 5,(0, r), {x*,y*) ^ ix**,y**) such that 

W'{y*,X) = W''{x*,X) and W'{y** , X) = Wix** , X). 

From (!25|) we have W'^{ym,X) = W'^{xm,X) = (xm,l/m,A) for m = 1,2. From point 4. 
in Theorem [TT] follows that 

Q,fo™ (r/, o W^iy*, X)-f],o W^iy**, A)) = 

Q^fo.w ((x*,i/*,r/,(A)) - (x**,i/**,r/,(A))) < 0, 

g^ba.k (77, o iy"(a;*. A) - r/, o A)) = 

Q^bac. iix\y\v^W) - ix**,y**,r],iX))) < 0. 

which implies that 

al\\x* -x**f + bl\\y* -y**f < 0, (26) 

a]^ \\x* - x**f + b^ \\y* - y**f < 0. (27) 

From ([21]) and ([27]) (keeping in mind that > 0, b^ < and that ajj < 0, feg > due 
to the reversion of the roles of x and y for the inverse map) follows that 

f II * **||2 ^ b II * **||2 ^ 7b II * **||2 . if II * **||2 

aoWx -X II >-ao||x -x \\ >bo\\y -y \\ > -bo\\y - y \\ , 

which contradicts (126]) . 

We now define xW '■= (^^cl/cA) for xq, yo such that W^{yQ,X) = W^{xo,X). By 
above arguments we know that x is a properly defined function. We need to show 
that this function is continuous. Let us take any A*, A** G ri^^{Bc{rji{\k), p)) for some 
{i, k) G J. From point 4 in Theorem [TT] follows that 

Q^fo™ (77, o x{X*) - fj, o xiX**)) < 0, (28) 
g^bac.(r/,ox(A*)-77,ox(A**)) <0. 
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Let us adopt notations fji o x(A*) = {x*,y*,9*) and fji o x(A**) = {x** , y** , 9**) . Note 
that from the construction of x follows that r]i{X*) = 9* and r]i{X**) = 9**. From f l28|) it 
follows that 

(4 + a^o) W^* - ^**f + W + &o) \\y* - y**f (29) 

< - (4 + Co) \\9*-9**\f 

= -(4 + ca lh.(A*)-r].(A**)ir 

From ([2lD it follows that + a}^ = \al\ - \a^^\ > and 6^ + 6|5 = - + > 0. By 
the fact that r]i is continuous and the fact that Cg < and Cq < 0, from (J29l) follows the 
continuity of x- 

We will now show that for any p G iy*(i?s(0, r) x A), (p) converges to x(^) ^-s n 
goes to infinity. Let us consider the limit set of the point p 

^{f^,p) = lim 12'' (p) = Q for some oo}. 

If we can show that u){ffp,p) is contained in (1 = x(A), then this will conclude 
our proof. We take any q = limfc_j.oo f^'^ip) from co{ffjj,p). By continuity of we know 
that q G W^. Suppose now that q ^ W^. This would mean that there exists an n > 
for which /^"(g) ^ 5„(0,r) x 5^(0, r) x A. Since 

hm /^-"(p) = 

we have that f^^{q) G uj{f^,p), but this contradicts the fact that uj{f^,p) C i?u(0,r) x 
5,(0, r) X A. 

Showing that all backward iterations of points in iy"(i?„(0,r) x A) converge to 
x(A) is analogous. □ 

Remark 14 . Let us note that during the course of the proof of Theorem [13] we have 
established more than just continuity of W^, and x- From our construction we 
know that for i G / 

V.oW''{x,v;\9))= {x,wnx,9),9), 
V,oW'{y,rj;\9)) = {wt{y,9),y,9) , 

v^oxiv^\o)) = Md),o), 

for continuous : i?„(0,r) x Be ^ 5,(0, r), wf : Bs{0,r) x Be ^ i?„(0,r) and 
Xi : B(. ^ Bu{0,r) x Bs{0,r). The inequality fl20|) from Theorem [TT] can be used to 
obtain explicit Lipschitz bounds for functions w^, wf. Also estimates ( 129|) can be used 
to obtain Lipschitz bounds for Xj. This means that we can get Lipschitz estimates for 
the invariant manifold x(^) together with Lipschitz estimates for its stable and unstable 
manifold. 
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4. Verification of covering and cone conditions 

In this section we show how covering relations and cone conditions can be verified with 
the use of local bounds on derivatives. The idea is to develop a simple automatised 
scheme which could be apphed in computer assisted proofs. In our approach we set up 
our verification so that we do not need to compute images of large sets (which in case of 
rigorous numerics is always troublesome). The scheme is based on iterates of a number 
of single points, combined with estimates on derivatives around their neighbourhoods. 
For any set V" C M" we define the interval enclosure of the derivative of / on \^ as 

WiV)] := 



mf - — (x),sup- — [x] 

x&v dXj x&v dXj 



for all i,j 



^5 



Let Ui-^ , [/j2 C A be such that dom/jj is nonempty. Assume that for any {c + u + s) x 
{c + u + s) matrix 

A e [dfi,i,{domfi,i,)] (30) 

we have the following bounds 

sup : \\vj\\ = 1} < Aij (31) 

inf {Pij^ill : ll^^ill = 1} >Aij, 

with i, j G {1, 2, 3} and f i, f 2, v^, representing the variables x, y, 9 respectively (note that 
Aij,A^j depend on the choice of 22^1)- In this section we shall use the bounds ( 15T]) for 
verification of covering and cone conditions. 

4.I. Verifying covering conditions 
We define a 3 x 3 matrix Ti^i^ as 
Ti^ii := (^ij)ij=i,...,3 

tu = All, ti2 = —A12, ti3 = —^13, 

^21 = ^21, ^22 = ^22, ^23 = ^23, (32) 

^31 = ^31, ^32 = ^32, ^33 = ^33- 

We will use notations R = (r„, r^, Vc) G and for q = (x, y, 6) eMy- xW x W and 
write 

N{q,R) := A^(g, r„, r^, rj. 

fi i 

We give a lemma, which can be used in order to verify that A^i =0 A'"2. 

Lemma 15. Let e > be a small number. Let Ni = N{qi,Ri) C dom/jjij be a ch-set. 
IfforRi = (r^,r2,r2) := Ti^^^Ri + {-e, e, e) we haverl,rl,rl > and for q2 := fi^iMi) 

Ikxg2|| +rl<l, ||vryg2|| +rl<l, heq2\\ +rl<l, (33) 
then for N2 := A^(?2, R2) we have Ni =0 N2. 
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A := Dfi^i^{qi + t{q - qi))dt G [rf^i^ (dom^ij], 

^0 
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we have estimates 

Ikx(/i2ii(g) - 92) 



kx(/i2n(g) - fi2iMi))\\ 

T^x (^j^ Dfi^^, (gi + t{q - qi))dt ■ {q - qi)^ 
= hxA{q-qi)\\ 

= \\Aun,^{q - qi) + Ai2TTy{q - qi) + v4i37re(g - qi 
>A^y^-M2rl-Mzrl 



~> r 

^ ' u-} 

hence ([2]) holds. Analogous computations for q ^ Ni give 

hy{fi2ii{<l) - 92)11 = hyMQ - <ll)\\ < ^2\r\ + ^22^5 + ^23'^c < ^s, 

lke(/»2ii(g) - 92)11 = heA{q - qi)\\ < A^irl + A^^r] + A^srl < r 
which proves ([3]). Conditions (!33|) ensure that C B. 



□ 



Example 16. We return to our Example [TOl The ch-sets from the example follow from 
Lemma [9] as Ni = N{0,Ri) where Rq = (r, r) and -R;+i = Ti^^^i^Rl + (—£,£,£) with 
T,,,, =diag«\4+i). 

Remark 17. When the x coordinate is strongly expanding, for practical reasons it might 
be beneficial to set significantly smaller than iTiTi^i^Ri. In such case the covering 

fi i 

Ni =^ N2 will still take place, but N2 will be a smaller set. This might give better 
bounds for next iterations of the map / and also keep the later constructed within 
B. Without reducing in the case when x is expanding, it might turn out that the 
sets Ni blow up quickly. 



4-2. Verifying cone conditions 

Now we shall present some lemmas, which will show how one can obtain condition ([1] 
from bounds on derivatives (|3T|) . The aim is to present a simple mechanism in which 
successive 7; are constructed. 

Let C = (cjj)jj=i_,,,_3 be a 3 X 3 matrix with coefficients 

Cii = Al^ - J2kf^l ^uAik C12 = ZlLl ^21^2fc 

C2I = A^2 " Ylk^2 ^12^1fc C22 
C3I = ii?3 - Efc^3 ^13^1fc C32 

(note that C depends on the choice of i2,ii)- 
We start with a technical lemma 



Yl'k=l ^22^2fc 
X]fc=l ^23^2A: 



Cl3 = ELi ^31 a 



3k 



C23 
C33 



ELi A32A3. (34) 



ELi A33 a 



k=l ^33^3A; 
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Lemma 18. Let 7 = (a, b, c) G M'^ and let A : — > ]^«+s+c ^ matrix for which 
the bounds l[31\) hold. // a > 0, 6 < 0, c < then for any p = (^1,^2,^3) G M'' x M" x 

Q^{Ap) > Qc^ (p) . 

Proof. Using the estimate 

±2 {AkiPi,AkjPj) > -AkiAkj {\\pi\f + \\Pjf) 

we obtain 

QMp) 

333 

= a ^ {AuPi, AijPj) + & 5Z (^2iPi, ^2jPi) + c (^3iPi, ^■ijPj) 

i,j=l i,j=l i,j=l 



= aY.MuPiW'' + bJ2M2^P^\\^ + cJ^MsiP^W^ 

1=1 i=l 1=1 

+2 ^ a {AuPi, AijPj) + 2^6 (^2iPi, ^2iPj) + 2 ^ c {AsiPi, Asjpj) 

i<j i<j i<j 

> \\Pif {aAl^ + bXi + cAIi) + \\P2f (a^?2 + + c^L) 
+ \\p3f{aAl, + bAl3 + cAl3) 

-aY^A.A, {M\' + \\p,\\') + bJ2W2, {M\' + \\p,\n 

i<j i<j 

+cJ2^3iA-ij (|b^||'+ Ibjil') 
i<j 

= iC^)l\\Plf+iC^)2\\P2f + {C^)Jp2f. 

□ 

Now we give a lemma which will be the main tool in the construction of 7; from 
Definition |9l 

Lemma 19. Let Ui^^Ui.^ C A and let N be a ch-set N C dom(/j2jJ. Let e > be a 
small number. Let C be defined by and e > 0. Assume that C is invertible and 
define 

Gi2h = C ^. (35) 
If for y = (a, b,c) := Gi^ij^'y+{e, e, e), we have a > 0, and b,c < then for any qi,q2 G N 
QyihhiQi) - fi2ii{<l2)) > Qj{qi - q2)- 

Proof. For 

^ := / Dfi2h{(l2 + t{qi - q2))dt G [rf/j^ij (dom^i J] 
^0 
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applying Lemma [TS] gives 

<5v(/i2n(?l) - 4*1 fe)) > QG,2n7(/i2ii(?l) - /i2iife)) 

= Q-yiQl - Q2)- 

□ 

Example 20. We return to Example [TOl The cones 7; follow from Lemma [19] as 
70 = (1,-1) and7;+i = {1 + e, (l + e)'^) -Gi^^^iji with d^^^i^ = diag ( (^m)^ ' (^wy) ' 
where ■ stands for the scalar product. 

4-3. Setting up local maps 

In this section we shall introduce conditions, which would ensure that the assumptions 
from Section 13.21 hold. Below we give a Lemma which will ensure conditions (|T^ and 

m- 

Let us note that in some cases conditions f lT^ and f|T5l) will follow from easier 
arguments or directly from the setup of the problem. Such is the case in our example 
from Section |5l 



Lemma 21. Let m > 1, A > and p > y + A. Assume that 

(i) for any l & I and any A G f/^ there exists a such that {t, /t) G J and 

MX)-7],{K)\\<A, (36) 
(a) for any 9 E Be and any i E I there exists an l E I such that 

{e, y^^) ^ ^ ° ^r') , (37) 

7],or]r\e)eB,{{],R-p-A). (38) 

For C^i defined as in ^^3^ , constructed for [d{j], o f]~^){(iom.{ri^ o ?7~^))] we assume 
that it is invertible and also that for 7 = (a, b, c) = C^^ji we have 

ao > ma, bo > mb, Cq > mc. (39) 

// assumptions [H [3 hold, then for any horizontal disc h in a ch-set with cones 
(M,j^) and for any i E I there exists (t, k) G J such that h(_B„(0,r)) C dom{fj^ o r/~^). 
Also for any qi, q2 in dovnifj^ofjl^) such that Q^_^{qi—q2) > we have p^ . Furthermore 
condition [TEi) holds. 



Proof. Let h be a horizontal disc in a ch-set with cones (M, 7^^). Take 60 = 7r0(h(O)). 
For any x G Bu{0,r) we have Q^-^{h{x) — h(0)) > 0, which implies that 

air^ > ai ||7r,,(h(x) - h(0))f > -Ci WMH^)) " ^of , 



Computer Assisted Proof for Normally Hyperbolic Invariant Manifolds 19 

hence 7re(h(i?u(0, r))) C Bc{9o, y^'^^r) fl B^- Taking l from assumption 2. for 6 = Oq, 
condition ( ^7]) implies that h{Bu{0,r)) C dom(?7t o ?^r^) and also 

\\v.°v7\0o)\\<R-p-^. (40) 

Take now any qi,q2 in dom(r7i o f^r^) such that Q-y^{qi — ^2) > 0. Applying ( l39l) and 
Lemma [19] gives 

<57o(^^ ° - ° ^^(^2)) > mg^(?7, o fil^iqi) -fj.o fil\q2)) 

= niQ^^iqi - q2) (41) 
>0, 

which proves f|T^ . Applying the bound in fHT]) for gi = h(a;i), g2 = h(x2) gives 

g^„(h'(xi)-h'(a;2)) >0, (42) 

which means that to prove ( |T5l) it is sufficient to show that h'(5„(0, r)) C M^^k for some 
K. Let A = r]~^{9o). We now take «; from assumption 1. For any x e S„(0,r), by (H2!) 
we have 

aor2>ao||7r,(h'(x)-h'(0))f 
> -Co||7r,(h'(x)-h'(0))f 
= -Co||7re(h'(x))-77,(A)f . 

This means that 



noih\B^iO,r)))cBMW,rJ^) 

V -Co 

hence 

\\7ie{h\x))-r],{X^)\\ < \\7re{h'{x))-r],i\)\\ + \\ri,{\)-V.{K] 



-Co 

< P, 

which gives h.'{Bu{0,r)) C Mj^^- What needs to be verified last is whether M^^k C B. 
From our construction vreMt^^ = Bc{rj^{X^), p). For 6* G Bc{ri^{X,^), p), using ( l36i) and 

ll^^ll < ||^^-7;.(A.)|| + ||r^,(A.)-r/,(A)|| + ||r;,(A)|| 

= \\0-v.iK)\\ + \h{K)~v.m + lh.°^r'(^o)|| 

<p + A + (i?-p- A), 
hence hqM^ G B^. □ 
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4.4- Normally hyperbolic manifolds from bounds on derivatives 

In Section 14.11 we have shown how covering relations from the chain (fT6l) can be 
constructed using bounds on derivatives of local maps. In Section 14.21 we have shown 
how the cones can be set up, using bounds on derivatives of local maps, so that the 
condition (ITTl) holds. Here we shall combine these results together in Theorem [131 

We shall use the notations T^jj^ and Gi^i-^ introduced in Sections 14. H 14.21 through 
equations ( 132|) . ( 134]) and ( 135|) . We will also assume that the assumptions from Section 
l3.2l hold. Here we introduce a definition which contains conditions which can be verified 
using computer assistance. We will later show that the conditions imply cone conditions. 

Definition 22. Assume that for any {lq, kq) G J there exists an n & N, a sequence 
Lq = Zq, ii, . . . ,in = Li and Hi such that {ti, Ki) G J and for 

= (x-,2/™,n:=/w.-.°---°/n.o(0,0,r/,„(A,,)), 

7 .= [a ,b ,c ) ■= Cri^i^^^ o ■ ■ ■ o Gt.iolo 
with m < n we have 

C + < 1, c + hni < 1, C + ll^"^ll < 1, 

<>r+||a;"||, r," + < r, (43) 

and 

a"^ > 0, > 6"^, > c'", 

a"- > ai, 6" > bi, c" > Ci. 

Then we say that f satisfies forward bounds. 

Remark 23. To verify that / satisfies forward bounds on needs to compute g*^, i?™" and 
7™. Let us note that in the case of it is enough to obtain bounds on a finite number 
of successive iterates of a single point. We therefore do not need to obtain bounds on 
images of large sets, which in practise would accumulate large errors. The and 
7"^ are constructed using local bounds on derivatives and are easily computable with 
computer assistance. Let us also note that to verify forward bounds we do not need to 
compute the composition function /" or its derivative (this would most likely cause big 
difficulties for high n due to complexity of such computations and also due to the fact 
that errors would accumulate quickly). 

Lemma 24. // / satisfies forward bounds then f satisfies cone conditions. 

Proof. We take any {lq, kq) G J, a sequence Lq = io,ii, . . . ,in = ti and an index /ti such 
that {li,Ki) G J from Definition [22l We define = := {r,r,p) and 

-R™ '.= Tj j -.RT ^ ~\~ { — 

iV^:=iV(g™,/?™). 

By (H3l) . taking e sufficiently small, we will ensure that C B. By Lemma [15] we 
obtain N^-i =^ ^ for m = 1, . . . , n and Nn ==^ M. 
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Figure 6. Misleading numerical plot of the attractor for T, obtained using double 
precision (grey), and the true invariant curve computed with 128bit accuracy (black). 

Now we define 7° = 7° := 7o and 

Taking e > small enough and applying Lemma [T^ we obtain f lT7|) . □ 

From now on let us assume that / satisfies forward bounds with 7q = 7o°™. 

Definition 25. Let -yl'^^ = (a{j, b|j, c|j) G be such that aj^, cjj < and > 0. l^e 

say that f satisfies backward bounds if f~^ satisfies forward bounds, with reversed roles 
of the X and y coordinates. 

Theorem 26. Assume that f satisfies forward bounds for 7o°™ = (aQ,bQ,CQ) and 
backward bounds for j^^'^^ = ^slq, bg, Cq) . // in addition inequality holds then there 
exists a normally hyperbolic invariant manifold inU, together with its stable and unstable 
manifolds W^, W^. 

Proof. This follows directly from Lamma [2l] and Theorem [131 d 

5. Example of applications 

Consider a driven logistic map 

T : 5^ X M ^ 5^ X R, 

T{e,x) = {e + a,l-a{e)x^), a(^) = ao + £sin(27re) (44) 

which differs from the well-known logistic map in the fact that the parameter a has 
been replaced by oq + e sin(27r6') and 6 has a quasiperiodic dynamics. Concretely we 
consider the parameter values oq = 1.31, e = 0.3 and « = 200' ^^ere g is the golden 
mean g = hence the dynamics on the base of the skew-product is slow. Numerical 

simulations in double precision (say, with mantissa of 52 binary digits) suggest that the 
map possesses a chaotic global attractor (see Figure [6l grey) . We will prove that this 
guess is not correct. When the same simulations are done with multiple precision, one 
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can guess that the attractor consists of two invariant curves (see Figure El black). We 
will use the method introduced in the previous sections to prove that T possesses a 
contracting invariant manifold and, in particular, that the red plot from Figure [6] do not 
shows the true dynamics. The same example was considered for other values of a and 
in a non-rigorous way in [2] to illustrate that one has to be careful with the arithmetics 
in simulations. 

5.1. Explaining the observed behavior 

To explain the reasons of the observed behavior it is worth to mention that in the 
example the parameter a of the logistic map ranges in [oq — e,aQ + e] = [1.01, 1.61]. For 
that range the attractor starts as a recently created (at a = 1) period-2 sink, followed by 
the full period-doubling cascade. Then one finds from several-pieces strange attractors 
to a single piece, interrupted by some periodic sinks and its corresponding cascades. 
When a moves with 6 one can question which is the "averaged" behavior. In particular 
the period-2 orbit is only attracting until a = 5/4. 

0.5 


-0.5 



-1.5 

0.2 0.4 0.6 0.8 1 




Figure 7. The integrand h{9) in for the parameter values: ao — 1.31, e = 0.30. 

To this end we can consider what happens for "frozen" values of a, denoting as Ta 
the corresponding logistic map. The orbit of period two, xi{a),X2{a) is given by the 
solutions of x^ — x/a + (1 — a)/a^ . In particular 

x^{a) = (1 - V4a-3)/(2a). (45) 

The differential of on it is 4(1 — a). To average with respect to 9 along the range, 
and noting that a — 1 > for the full range, we have to consider the average of the 
Lyapunov exponent given as 

\ [ log(4(ao-l + £sin(27r^)))d^^ = ^log(2(ao-l + v/(ao-l)2-£2))(46) 

which for ao = 1.31,5 = 0.3 gives Aqo ~ —0.12666931. The integrand is shown in Figure 
[71 For the skew product, assuming a ^ Q and sufficiently small the two curves which 
form the attractor, as will be proved later, are very close to the curves Xi{a),X2{a) of 
the frozen system. Figure El displays the lower one. Also the Lyapunov exponent of the 
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driven map with a = g/N, N = 200, computed using 10^ iterates after a transient also 
of 10^ iterates is A200 ~ -0.12680. Using other values of A^, like 100,400,800, 1600 the 
respective values A^v obtained are -0.12725, -0.12670, -0.126696, -0.126689, tending 
to the limit Aqq. 



-0.1 

-0.2 

X 

a 

U 0.2 0.4 0.6 0.8 1 



Figure 8. The lower part of the attractor, the graph of xi{a{9)), for the parameter 
values: gq ~ 1.31, e ~ 0.30. 

The numerical difficulties are easy to understand. To compute the Lyapunov 
exponents, starting at a point xq and an initial vector vq = 1 and setting 5*0 = 
we compute recurrently 

%+l = DTa{Xj){Vj), Xj+i = Ta{Xj), Uj+i = \Vj+i\, Vj+i = Vj+i/Uj+i, 

Sj+i = Sj + log(nj+i). 

The values Sj are denoted as Lyapunov sums and the average slope as a function of j 
(if it exists) gives the Lyapunov exponent A. For details and generalisations see, e.g., 
pT] and [H] and references therein. 

Even when A is negative it can happen that partial sums have strong oscillations. 
Given the values of Sj,j = 0, . . . ,k let (5*^)111111 be the minimum of these values and 
introduce Ok = Sk — (5'fc)iiiin- We define the maximal oscillation of the Lyapunov 
sums as OS = max{Ok}. The Figure [9] shows the behavior of Sj for a = g/200 
and also some of the initial oscillations for a = (7/I6OO. A non-rigorous computation 
of OS for = 100,200,400,800,1600 with 10^ iterates after a transient gives the 
values 28.845, 56.761, 112.632, 224.379, 447.874, respectively. This implies a loss in the 
number of decimal digits equal to these values divided by log(lO). In particular, between 
24 and 25 digits for A^ = 200, which explains the failure seen in Figure [61 For small a 
the maximal oscillation tends to be 
1 f'' 

- / h{9)d9, (47) 

where h{6) is the function which appears as integrand in fHB]) and it is extended by 
periodicity outside [0,1] while ^1 = | - ^cos~^(0.2), 62 = I + ^cos"^(0.2) are the 
values at which h becomes equal to zero (see Figure [7]). The value of the maximal 
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oscillation in (gT]) is ^ 0.172660185/a for small a, that ifa = g/N becomes ^ 0.27937A^ 
in good agreement with the previous results. 
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Figure 9. Oscillations of the Lyapunov sums. Left: the Lyapunov sums for N = 200. 
Right: some initial sums for N = 1600. Parameter values: oq — 1.31, e = 0.30 and 
a = g/N. 

Using these ideas one can even predict when we shall observe that the attractor 
produced by simulations with not enough digits seems to indicate that it is not a period- 
2 curve. Assume that we do computations with d decimal digits and that in a plot like 
the one in Figure [6] one can distinguish pixels which are a a distance of 10"^. In our 
example reasonable values of d,p are 16 and 4. This means that from 62 — 1, when h 
becomes positive, till some unknown 6^ when the "departure" of the iterates from the 
curve become visible, the factor of amplification of errors is 10*^"^ or, in logarithmic 
scale {d — p) log(lO). This requires 
1 f'^ 

- / h{e)de = (rf-p)iog(io). 

In our example one finds 9d ~ 0.258 in good agreement with the observed numerics in 
Figure |6] . In a similar way one can predict the "landing" value 9i at which the points 
seen as chaotic in Figure E] are close enough to the real invariant curves. As the distance 
from the chaotic points to the true attractor is of the order of 1, the condition is now 
1 f'^ 

For the example one obtains 9i ~ 0.629, again in good agreement with the observed 
numerics. 

This "delayed" observation of the expanding and compressing regimes is similar, 
but now due to purely numerical reasons, to the delay of bifurcation that can be observed 
in systems depending on a parameter which has slow dynamics (see [16] and references 
therein) . 



5.2. Some limit cases 

Now we discuss two limit cases. First one is the case in which a{9) covers a wide range. 
Second one aims at describing the differences between the union of the curves Xi{a{6)) 
and X2{a{9)) and the true attractor for a small enough. 
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According to ( I46p and assuming that for a sufficiently small the attractor is close 
to the union of the curves Xi^2{ci{d)) it is enough to take Oq = 1.5 — 6i, e = 0.5 — 61—62 
with < 62 < 6f to have a negative limit averaged Lyapunov exponent Aqo. If 61 is small 
the values of a almost cover the full range (1,2). The Figure [TO] displays results of the 
observed behavior using double precision for the values 61 = 0.005, 62 = 10~®, a = 
^/60000. 
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Figure 10. Simulations in double precision for values of oq, e such that a{9) almost 
covers the range (1, 2) and a very small. See the text for the numerical values used. 

The figure is reminiscent of the "bifurcation diagram" of the logistic map. In fact, 
a typical way to compute the diagram consists of taking a sample of values of a, do some 
transient iterates and display some of the next iterates. Now the value of a is changed 
at every step according to dH]) but very slowly, and the transient is discarded. From 
6' = 3/4 (for which the minimum of a{9) is achieved) to 9 = 5/4 (mod 1) (for which 
one achieves the maximum) the plot looks like that diagram, except for the bifurcation 
delays at the period doublings from period 2 to period 4 and successive ones. In the 
range 6 e [1/4, 3/4] the reverse situation is seen, but now with much smaller bifurcation 
delays. The authors do not know if, for the present values of the parameter, the attractor 
will become close to the union of Xi{a{9)) and X2{a{9)) for computations done with a 
huge number of digits. 

To look for the expression of the attractor as the union of two smooth curves, 
assuming it is of that type, we restrict our attention to the lower part of it, close to 
Xi{a{9)) as given in (l45l) . In principle it is convenient to work with but, as the 
eigenvalues of along the points of period are negative, we prefer to work with T^. 
We look for the attractor as the graph of a function expanded in powers of a 

G{9) = Go{9) + aGi{9) + a'G2{9) + . . . , (48) 

where Gq{9) = xi{a{9)) is the zeroth order approximation. The map T^{9,G{9)) is 
0{a) close to the identity. Hence, it can be approximated by a smooth flow (see [I] for 
proofs, an example of application and additional references, as well as [15] for general 
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results) and the curve we are looking for is a periodic solution of this flow. But we shall 
proceed by imposing directly the invariance condition. 

Starting at a point of the form {9, G{9)) and doing four iterations using the values 
a{6), a{6 + a), a{0 + 2q;), a{6 + 3a) we should have 

T\e, G{e)) -{6 + 4a, G{e + 4a)) = 0. (49) 

Given values of ao, e it is a cumbersome but elementary task to obtain in a recurrent 
way the expressions of Gi, G2, ■ ■ ■ from pOl) . It is essential to reduce the dependence 
in Gq{6) using the equation satisfied by xi(a) to decrease the order of the powers of 
Go which appear to just the first one. We note also that in the computation of all the 
terms Gj there appears 16a^ — 32a + 15 = (4a — 5) (4a — 3) in the denominator, which 
cancels for a = 5/4, but a careful examination allows to show that the factor 4a — 5 is 
also present in the numerator. 
In this way one obtains 

3 - 2a - (8a - 9)/V4a - 3, 



^'^^^ 2a2(4a-3) 
where a stands for a{9) as introduced in (l44l) . 



-2TTe cos{2tt6), 



(50) 
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Figure 11. Graphs of Gi{0) (left) and G2{e) (right) for oq = 1.31, e = 0.30. 

The computation of G2 is much more involved. The simplest expression is given 
as a rational function depending on a{9), Gq{9), Gi{6) and up to the second derivatives 
of these functions with respect to 9. Instead, Figure [TT] displays the graph of Gi and 
G2 for ao = 1.31, e = 0.30. The graph of Go{e) is very close to the attractor shown in 
Figure [H 

To see tiny details on the attractor Figure [12] displays the differences between the 
lower part of the attractor, computed with enough digits, and the approximation in (HHj) 
up to order 2 in a. 

The left part shows tiny oscillations which were not visible in Figure [81 They reach 
a maximum at the value 6 = 61 for which h{6) in (H71) changes from positive to negative. 
As one can expect the shape of these oscillations is a bump function multiplied by a 
periodic function (close to a sinus) with period 4a. A similar behavior is observed for 
many other values of aQ,e and a. When the oscillations start at a larger distance from 
61 they can amplify is such a way that the attractor is no longer the union of the two 
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Figure 12. Differences between tlie attractor and the second order approximation for 
flo = 1.31, e = 0.30 and a = g/N. Left: N = 200. Right: N = 1000. 

curves. One can suspect that it becomes a non-chaotic strange attractor (see, e.g., [10] 
and [12]). In contrast, with the same values of ao,e but for = 1000 the oscillations 
are not observed and the very small differences in the plot on the right hand side of 
Figure [12] are mainly due to the third order term in ( HHl) . 

5. 3. Computer assisted proof of existence of invariant curves 

In this section we apply our method from Sections 131 15.3. II to prove that for parameters 
= 1.31, £ = 0.3 and a = with g = ^"-^ the map T has an invariant curve. 
Around a neighborhood of the numerical guess for the attractor, the map is 
locally invertible. This is due to the fact that our curve is separated from the x-axis. 
For our proof we consider 

/ = T-\ 

The attractor is first computed (nonrigorously) by iterating T forwards in time. We 
then choose a set V arround the attractor (see Figure [131 gray). For most 9 the set is a 
0.001 radius neighbourhood of the attractor. Close to the angle 6* = | we choose V to 
be tighter, so that we are sure that it lies within the domain of / (see Figure [T3l) . Our 
aim is to prove that inside of V we have an invariant normally hyperbolic curve of /. 
The map / is not uniformly expanding in the x direction. Over one part of the set V 
the map / is strongly expanding, elswehere it is contracting. A part of the expansion 
region, which we denote as W C V, is depicded in red and green (the green region is on 
the left tip of the red region and is poining towards the attractor) on Figure [13] On this 
set we place ch-sets Ni, . . . , Niqs of width |, starting with Ni on the left and finishing 
with A^i68 on the right. We shall use a notation 

I 

i=k 

Our ch-sets are parallelograms. The coordinate x is globally expanding for / and 
coordinate is normal (our map does not posses a globally contracting coordinate y). 
The exits sets for the ch-sets are the top and bottom edges of the parallelograms. 
The map / moves the ch-sets to the left. We distinguish two parts of the set IA\ the set 
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Figure 13. Positioning of our ch-sets (green and red) and the set V (gray) relative to 
the attractor (on this plot on the f?-axis in black). 

I ' ' 1 0.002 I ' ' ' 1 




Figure 14. The ch-sets: N-^ , ■•^4 in green and ,...,N^Qg^ in red (plotted 
relative to the attractor), together with f{N^),...,f{N^Qg) in dark blue and 
/i28(7Vr ),..., /i28(iV4-) in black. 



Ui^4 in our plots is denoted in green colour, U^^ies is denoted in red. Since the width of 
the ch-sets is f, for A; G 5, ... , 168 we have 

T^efiNk) C neNk-i. 

In Section 15.3.11 we shall show that (see Figures [HI [T5l) 



Nk^Nk^4 for A; G {5, . . . , 168}, (51) 



and also that for i = 1, 4 (see Figure [T] 

fl28 

In Section 15.3.21 we show how to verify cone conditions. In Section 15.3.31 we briefly 
discribe the tools that were used to conduct the proof. 

5.3.1. Verification of covering conditions To describe how covering conditions are 
verified we start with a seemingly unrelated discussion. Consider a polynomial p : 
[0, r] — )■ M of degree n 

n 

p{e) = Y,a,e^, 

j=0 
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Figure 15. Closeup of the covering Ni Ni^4 for i = 5, 6, 7, 8 (plotted relative to 
the attractor). 



and a function : R — )■ R. Using Taylor expansion and defining two polynomials p and 
p, of degree n 

n-l 



m-9opio) + j:{^'-^io))e^ 



(52) 



1 / d^igop) 1 cP+M^ 

H r 77 H sup -— — (v)w , 

n\\ d9" ^ ^ n + 1 .^^J,,] d9^+^ v ; | , 

n-l 

/ I i^./ I n Ti 1 

(0) 0' (53) 



for any 6 E [0, r] we have 

m<g{p{9))<m- (54) 

For any i = 1, 168, the exit set consists of two lines and can be expressed 
using two polynomials (in fact these are affine functions) Pi,pf : [0, |] — ?► R, pf{0) = 
afo + af^^O, Pi{6) = a"o + a^iO and a point G [0, 1), 

AT- =iV-UiV-, 

pf{9)<pn9) for^^G[0,|]. 
We will now show how to construct a ch-set M such that 

Ni ^ M. (55) 
We first verify that for any point {6, x) G Ni we have f^(x, 6) < 0. We then take 

g'^iO) := /(g. + ^,pfW), := f{q^ + e,pm), (56) 



Computer Assisted Proof for Normally Hyperbolic Invariant Manifolds 



30 



and construct p^{9) = p{0) using (13^ and p'^iO) = p{0) using (15^ . taking g as functions 
gf" and 5''^ respectively. Formula fl54|) guarantees that f{N~^) lies above the graph of 
p"(6') and that f{N~^) lies below the graph of p'^{6). If we now set 

M- = M-UM-, 

M,- = {(/(^), - 2a + G [0, |]}, (57) 
M- = {(p"(e),g.-2a + ^^)|eG[0,|]}, 

and take M to be the set of points which lie above and below M~ then ( 155|) holds. 

For i = 5, 168, after applying the above procedure to abtain M which is covered 
by Ni, we compute bounds on the images of sets 

forj = 0,...,9, (58) 

in local coordinates of ch-sets iVj_4, to verify that we have flSTl) (subdividing [0, |] into 
ten intervals turns out to be sufficient for all i G {5, 168}). 

For i = 1,...,4 we need to iterate the procedure (l57Il many times to obtain a 
sequence of covering relations 

AT, ^ Ml ^ M2 ^ . . . ^ Mi27 ^ t/5,168. 
During our construction we make sure that all sets Mk for k G {1, 127} lie in V, which 
readily holds since the sets are very strongly contracted. Each covering Mk =^ M^+i 
holds by construction. Verifying that M127 =^ ^5,168 is done analogously to (135]) . 

In our computer assisted proof we take the degrees of polynomials for the edges of 
the sets Mk as nine, which means that we need to perform computations. Let us note 
that computationally this is not as heavy as might seem, since the C^^ computations 
are performed for one dimensional functions g^{9) and g'^{9) (see (1561) ). The reduction 
of dimension truly pays off, since the difference between computations in one and 
two dimensions is substantial. 

The estimates obtained by us are very accurate. In Figure [I6] we give a plot of 
M^28uy which is the lower bound estimate of the image of N^^ after the final step 
in our procedure (in black), and compare it with ten points from N^^, iterated non- 
rigorously with high precision computations (in red). The curve lies below the points, 
as should, but this is impossible to distinguish from the graph. The right hand side of 
Figure [16] gives the plot of the difference of the rigorous lower bound and non-rigorous 
computation. They turn out to be very close. 

Remark 27. The high order computations and multi-precision in current approach seem 
essential. The sets Mk constructed with our procedure are very strongly contracted. The 
distance between the two curves of Mj7 at the tightest spot is of order 1.125 x 10~^^, 
which is extremely thin when compared to the width of the curves ^ ~ 1.545 x 10"'^; 
and yet, with our approach, with little effort we are able to rigorously keep them 
apart. Any standard approach, such as performing C° computations on sets or careful 
linearization with techniques through local coordinates, is likely to fail. 
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Figure 16. Left: Rigorous bound on the image of an edge of one of ch-sets after 128th 
iterate of the map (in black), together with non-rigorous computations using multi- 
precision (red). Right: The difference between rigorous lower bound and non-rigorous 
computations. 



Remark 28. We believe that using a "parallel shooting" type approach it should be 
possible to conduct the proof using double precision and computations only (for this 
we would need an good apriori guess for the position of the curve). Such approach 
could produce a rigorous-computer-assisted proof using double precision of an invariant 
a curve, which is not detectable numerically with double computations. This shall be a 
subject of forthcoming work. 



5.3.2. Verification of cone conditions To verify cone conditions let us first rescale our 
coordinates by 

Taking /3 sufficiently large, choosing sufficiently many points Aj G [0, (3) and taking 
hi := i (c« (Xi) - c'^ (Ai)) , := (A^, c'^ (A^) + and V := V n ([A^ - h,, A, + h,] x M) 
we can construct local maps 

for which fii{Vi fl c") = B^x {!}, fjiiVi fl c*^) = B^x {—1} and which are arbitrarily close 
to a linear map q — > -^{q — qi). In these local coordinates, by taking sufficiently large 
/3, we have the following bound on derivatives of local maps (assuming that we choose 
i,i and p such that p G dom{fij) ^ 0) 

Dfij = (?7i o 7^ o / o 7^ 1 o 777!) (p) 
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which in turn is arbitrarily close to -j^diag^-^, -j^). This means that by using the artificial 
rescaling 7^ (without the actual need to apply it in practice for our computer assisted 
proof), we can divide the region V into a finite number of sets Ui, . . . ,Un (V C Uili ^j); 
and verify cone conditions using interval matrices diag{ [^{Ui)] , [^(f^i)] ) and applying 
Lemma fT9l For our proof we take 70 = (a, b) = (1, —1) , which means that the quadratic 
form for our cones is simply 

Q.,X0,x)=x^-9^. 

If we take 71 = ((1 — £:), 1) for any small parameter e > then by choosing sufficiently 
large /3 Assumption [7] is satisfied (since any switch to new coordinates is arbitrarily close 
to identity). This means that we can take 71 = 70, provided that all the inequalities in 
our verification of cone conditions in the computer assisted proof are strict. 

5.3.3. Tools used for the proof Our proof has been conducted with the use of the 
CAPD library (http:/ /capd. ii.uj.edu. pi) developed by the Computer Assisted Proofs in 
Dynamics group. We have used the multi-precision version of the library running at 
128 mantisa bits accuracy (which is approximately equivelent to tracking 40 digits). 
The computations have been performed with assistance of the Flexible Automatic 
Differentiation Package FADBAD++ (www.fadbad.com). The proof takes 16 seconds 
running on a 2.53 GHz laptop with 4GB of RAM. 

6. Final comments 

In this paper we have presented a version of a normally hyperbolic invariant manifold 
theorem, which can be applied for rigorous-computer-assisted proofs. We have 
successfully applied our method to an example in which standard double precision 
simulations brake down and produce false results. This demonstrates the strength 
of our method, that it can handle numerically difficult cases. It needs to be noted 
that to apply our method we have used multiple precision for our computer assisted 
computations. For our proof we also needed to apply a high order method which relied 
on computations. We believe that it should be possible to devise a similar in spirit 
method, which would give proofs without multiple precision and using computations 
only. This will be the subject of our future work. 
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